###### 4.2 - Cervejas

cervejas <-c("chope","lata","garrafa","chope","garrafa", "garrafa","lata","lata","nenhuma","lata","garrafa","garrafa", "garrafa","lata","lata","lata","garrafa","lata","chope","nenhuma", "garrafa","garrafa","garrafa","chope","garrafa","garrafa","chope","garrafa","lata","lata")
freq.cerv <- table(cervejas)
freq.cerv
barplot(freq.cerv)

as.numeric(freq.cerv)

dotchart(freq.cerv, xlab= "Tipos de cervejas preferidas")

#o dotchart exibe maior informao com menos tinta, portanto apresenta maior razo dado/tinta.

###### 4.3 Caixetais

## Construa um histograma do dap1) dos fustes dos caixetais.
## Construa histogramas da altura das rvores para os diferentes caixetais ('local').

caixeta = read.table("caixeta.txt", header=T, sep=",")

head(caixeta)
dap <- caixeta$cap/pi
hist(dap)

levels(caixeta$local)
t.tot <- table(caixeta$h, caixeta$local)
t.tot
t.chauas <- t.tot[,1]
t.chauas
t.jureia <- t.tot[,2]
t.retiro <- t.tot[,3]

par(mfrow=c(1,3))
hist(t.chauas)
hist(t.jureia)
hist(t.retiro)
par(mfrow=c(1,1))


######## 4.4 Eucaliptos

####Neste exerccio, use o conjunto de dados Inventrio em Florestas Plantadas de Eucalyptus grandis.

getwd()
egrand <- read.table("egrand.txt", header=T, sep=";")

## 4.4.1 Utilize o grfico boxplot para analisar o DAP de rvores de E. grandis em funo das variveis regio (regiao) e rotao (rotacao).

#por regiao
levels(egrand$regiao)
dap.Bofete <- egrand[egrand$regiao=="Bofete",]
dap.Botucatu <- egrand[egrand$regiao=="Botucatu",]
dap.Itatinga <- egrand[egrand$regiao=="Itatinga",]
dap.Salto <- egrand[egrand$regiao=="Salto",]

dap.reg <- list(a=dap.Bofete$dap, b=dap.Botucatu$dap, c=dap.Itatinga$dap, d=dap.Salto$dap)

boxplot(dap.reg, names=levels(egrand$regiao))

#por rotacao
str(egrand$rotacao)
range(egrand$rotacao)
summary(egrand$rotacao)
unique(egrand$rotacao)
dap.rot1 <- egrand[egrand$rotacao==1,]
dap.rot2 <- egrand[egrand$rotacao==2,]
dap.rot <- list(dap.rot1, dap.rot2)

boxplot(egrand$dap ~ egrand$rotacao, data=egrand)

boxplot(egrand$dap ~ egrand$regiao, data=egrand)

## 4.4.2 Avalie a normalidade da altura do conjunto total de rvores com um grfico quantil-quantil contra a distribuio normal.
hist(egrand$ht)
qqnorm(egrand$ht)
qqline(egrand$ht)


######### 4.5 Mais Caixetais 

caixeta = read.table("caixeta.txt", header=T, sep=",")
str(caixeta)
unique(caixeta$especie)

##clculo do dap
caixeta$dap <- caixeta$cap/pi

##selecionando apenas a especie Tabebuia cassinoides
Tab <- caixeta[caixeta$especie=="Tabebuia cassinoides",]

## tentativas de calcular a somatria dos daps (o certo seria calcular a rea de cada fuste, obter a rea total e depois encontrar o dap equivalente a rea basal)
arv <- which(Tab$arvore==Tab$arvore)
arv.dap <- tapply(Tab, arv, FUN= sum(Tab$dap/pi))
tapply(Tab, subset(Tab$arvore, Tab$arvore==Tab$arvore), FUN= rowsum(Tab[,-c(1,7)], Tab$dap))

## grficos
coplot(Tab$dap~Tab$h|Tab$local)
scatter.smooth(Tab$dap, Tab$h)

lattice::xyplot(Tab$dap~Tab$h|Tab$local)
